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Abstract 


Real-time quantitative PCR (qRT-PCR) is often used as an effective experimental method for analyzing 
gene expression. In this method, normalization of target gene expression levels must be performed using 
housekeeping genes (HKGs). HKGs are used to compensate for difference between samples due to diverse quality 
and quality of RNAs and different reverse transcription yield. For an ideal reference gene, constant expression levels 
across different samples of one experiment is necessary. In the current study, expression stability of four candidate 
references genes including Beta actin (ACTB), glyceraldeyde-3-phosphate dehydrogenase (GAPDH), hypoxanthine 
guanine phosphoribosyl transferase (HPRT1) and Beta-2-Microglobulin (82M) following retinoic acid (RA) 
treatment in embryonal carcinoma NCCIT cells were evaluated. NCCIT cells were exposed to RA (10 uM) for 14 
days to induce differentiation. RT-qPCR for candidate references genes was performed and normalization between 
untreated and RA-treated cells was performed using identical sample input amounts. Expression of OCT4, SOX2, 
NANOG during RA-induced differentiation was assessed by quantitative real-time PCR. RT-qPCR results indicated 
significant difference in expression level of GAPDH between untreated (Ct mean: 19.36667+ 0.28) and RA-treated 
(Ct mean: 28.94+ 0.18) NCCIT cells. However, transcriptional level of ACTB, HPRT and B2M remained unchanged 
after RA treatment. (RT-PCR analysis using ACTB, HPRT and B2M showed treatment of NCCIT cells with RA lead 
to significant down regulation of OCT4 (79%), NANOG (71%) and SOX2 (96%) transcript. ACTB, HPRT and B2M 
were recognized as valid reference genes for analysis of gene expression during RA-induced differentiation of NCCIT 


DOI: 10.22067/jemr.v12i2.88457 


cells, while GAPDH was not suitable. 
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Introduction 

Real-time reverse transcription quantitative 
PCR (qRT-PCR) has emerged as one of the most 
powerful tools for quantifying gene expression 
levels in a variety of gene expression studies. In real 
time PCR, it is required to normalize variations in 
quantity and quality of starting RNA and cDNA 
synthesis efficiency among different samples. For 
this purpose, internal control, often referred to as a 
housekeeping or reference gene are widely used 
(Nygard et al., 2007). Housekeeping genes 
(HKGs) are constitutive genes that are required for 
maintenance of basic cellular function (Richly et al., 
2003). Since, HKGs play a critical role in accurate 
normalization of qRT-PCR data, an ideal HKGs 
must be expressed ina _ stable and _ non- 
regulated constant level across different cell types 
and experimental situation (Chen et al., 2006; 
Turabelidze et al., 2010). 
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Beta actin (ACTB), glyceraldeyde-3-phosphate 
dehydrogenase (GAPDH), hypoxanthine guanine 
phosphoribosyl transferase (HPRT1) and Beta-2- 
Microglobulin (82M) are among the most commonly 
used HKGs in qRT-PCR studies (Dheda et al., 2004; 
Radoni€¢ et al., 2004). However, several reports have 
shown that these HKGs have variable expression 
levels in different experimental conditions (de Kok 
et al., 2005; Glare et al., 2002; Selvey et al., 2001). 
Therefore, there is no universal HKGs having stable 
expression in all tissues under all experimental 
conditions and it is clear that suitable HKGs must be 
selected for particular sample sets and experimental 
models. 

The purpose of our current investigation is to 
evaluate the expression stability of four candidate 
HKGs (ACTB, GAPDH, HPRT and B2M) during 
retinoic acid (RA)-induced differentiation of 
embryonal carcinoma (EC) NCCIT cells. According 
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to our results, GAPDH showed high variability of 
expression between untreated and RA-treated 
NCCIT cells. Therefore, although GAPDH is a 
common reference gene in many experiments, it 
cannot be an optimal reference gene in this 
experiment. On the contrary, expression levels of 
ACTB, HPRT and 82M remained unchanged in 
NCCIT cells after RA treatment, so they can be used 
as a suitable reference gene for analysis of gene 
expression in NCCIT cells exposed to RA. 


Materials and Methods 


NCCIT cell culture and RA treatment 

NCCIT cells (Pasteur Institute of Iran, Tehran) 
were cultured in RPMI-1640 medium (Gibco, Grand 
Island, USA) supplemented with 10% fetal bovine 
serum (FBS; Biowest, France) and 100 U/mL 
penicillin and streptomycin (Biowest, France) in the 
presence of a humidified atmosphere with 5% COz 
at 37 °C. According to many researches, to induce 
differentiation, NCCIT cells were treated with 10 
uM RA for 2 weeks (Jin et al., 2012; Kooistra et al., 
2009; Park et al., 2014). 


RNA extraction and cDNA synthesis 
RNA isolation kit (DENAzist Asia, Mashhad, 


Table 1. List of different PCR primers used in the study 
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Iran) was used to extract total RNA from untreated 
and RA-treated NCCIT cells. Afterward, quality and 
quantity of RNA were assessed via nanodrop and 
agarose gel electrophoresis. To avoid DNA 
contamination, total RNA was treated with DNase I 
(Cat. No. EN0521; Thermo Fisher Scientific). For 
cDNA synthesis, 1 ug RNA was reverse transcribed 
using 200 U of M-MuLV reverse transcriptase (Cat. 
No. EP0441; Thermo Scientific) in presence of 5 uM 
oligo (dt)18 (Cat. No. MANO013109; Thermo 
Scientific), 1 mM dNTPs (Cat. No. RO192; Thermo 
Scientific) in accordance with the manufacturer’s 
instructions. 


Quantitative Real-Time PCR 

Real time PCR was performed using real-time 
PCR system (Analytik Jena, Jena, Germany). For 
each reaction, 12.5 ul of SYBR Green, 0.5 ul of 
forward and reverse primers (10 pM) and 1p] cDNA 
in total volume of 25 yl was used. Real time program 
for all genes was: 95°C for 4 min followed by 40 
cycles of 95°C for 30 s, 62°C for 20s, and 72°C for 
10 s. To derive melting curves, at the end of the PCR 
run, temperature was increased in steps of 1°C for 10 
s from 61°C to 95°C. The sequence of primers and 
product length are described in Table 1. 


Gene nanie Sequence (5* to 3’) Product 
size (bp) 
POU class 5 homeobox | (POUSF1) F:CCGAAAGAGAAAGCGAACCAGTAT 145 
Alias symbols: OCT4 R: CCACACTCGGACCACATCCTTC 
Nanog homeobox (NANOG) F: AATACCTCAGCCTCCAGCAGATG 149 


R: CTGCGTCACACCATTGCTATTCT 


SRY-box transcription factor 2 (SOX2) 


F:GGGAAATGGGAGGGGTGCAAAAGAGG 151 


R: TTGCGTGAGTGTGGATGGGATTGGTG 


Hypoxanthine guanine phosphoribosyl 
transferase (HPRT) 


F: TTTGTTGTAGGATATGCCCTT 168 
R: ACATTGATAATTTTACTGGCGAT 


F: ACCACCTTCAACTCCATCATG 120 


Beta-actin (ACTB) 


Glyceraldehyde-3-phosphatedehydrogenase 
(GAPDH) 


Beta-2-Microglobulin (62M) 


R: CTCCTTCTGCATCCTGTCG 


F: GACCACTTTGTCAAGCTCATTTCC 151 
R: GTGAGGGTCTCTCTCTTCCTCTTGT 


F: CTCCGTGGCCTTAGCTGTG 69 


R: TTTGGAGTACGCTGGATAGCCT 
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Data analysis 


Fold change in expression of the HKGs in RA- 
treated versus the untreated samples 

Reverse transcription was performed with the 
same amount (1 pg) and optimal quality of total 
RNA from two samples (untreated and RA-treated 
NCCIT cells), so mean cycle threshold (Ct) values 
of two samples were used to evaluate expression 
stability of candidate HKG. 
In order to calculate fold change in expression of 
the HKGs in the treated versus the untreated 
samples, at first mean Ct values was calculated as 
2 and fold change is calculated as Mean 2M treated 
cell! Mean oS untreated cells. 


Fold change in expression of pluripotency 
factors during RA-induced differentiation of 
NCCIT cells 

To assess expression of OCT4, NANOG and 
SOX2 in RA-exposed cells in comparison with the 
untreated control cells, B2M, HPRT and ACTB 
were used as HKGs and data were analyzed 
by relative quantification using the comparative 
Ct method: Fold change = 2 ““*, 
AAC = [ACt treated cells (Meant Ct target gene -Mean 
Ct reference gene) _ ACt control cells (Meant Ct target gene 
a Mean Ct reference gene) | 
Statistical analysis 

In this experiment, NCCIT cells were cultured 
and treated with RA in three “6 cm cell culture 
plate” and then differentiated cells were merge and 
one RNA was extracted from them. The same was 
done for untreated NCCIT cells. For qPCR, there 
are three biological repeats for each gene. Mean Ct 
values are shown in Supplementary table 1. T-test 
was used tocompare expression of genes 
between treated and control NCCIT cells. P-value 
less than 0.01 was considered statistically 


Soltanian and Sheikhbahaei 


significant. 
Results and Discussion 


Identification of ideal HKGs that are stably 

expressed under various experimental conditions 
and tissues of interest would provide a powerful 
tool for normalization of target gene expression 
levels in (RT-PCR. On the contrary, poor selection 
of HKGs can invalidate the normalization process 
and lead to the generation of misleading 
information (Chen et al., 2006; Glare et al., 2002; 
Haller et al., 2004). Unfortunately, it has been 
proved that expression of commonly used HKGs is 
affected with treatment or physiological state 
(Bustin, 2000; Janovick-Guretzky et al., 2007; Wu 
and Rees, 2000). Therefore, studies to find suitable 
HKGs with a relatively stable expression level 
across many different cell type and conditions have 
received more attention in gene expression analysis 
(Vandesompele et al., 2002). In the present paper 
we have therefore looked at the expression stability 
of commonly used HKGs including ACTB, HPRT, 
B2M and GAPDH during RA- induced 
differentiation of NCCIT cells. 
For this purpose, normalization was performed 
between untreated and RA-treated samples using 
same amount of RNA (1 ug) in cDNA synthesis 
and optimal quality of RNAs were confirmed by 
nanodrop (260/230: 2-2.2 and 260/280: 1.8-2) and 
electrophoresis (supplementary Figure 1). After 
qRT-PCR, analysis of melting curves clearly 
indicated that each of the primer pairs described in 
Table 1 amplified a single expected product with a 
distinct Tm (supplementary Figure 2). The 
accuracy of the amplification reaction was 
validated by gel electrophoresis (Figure 1). 


Figure 1. The accuracy of qRT-PCR was further validated by gel electrophoresis. RT-PCR products from left 
to right are GAPDH: 151 bp; B2M: 69 bp; ACTB:145 b and HPRT: 168 bp genes for untreated and RA-treated 
NCCIT cells respectively. A molecular weight marker (100 bp ladder) is used. 
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Ct values of the candidate HKGs obtained by 
qRT-PCR in both samples (untreated and RA- 
treated NCCIT cells) were used to assess their 
expression levels and stability (supplementary 
Table 1). In comparison with other HKGs, the 
average CttSD (25.8+0.025) of B2M in untreated 
NCCIT cells was the highest, meaning it had the 
lowest expression level than other and the most 
abundantly expressed gene was ACTB and its 
average Ct+SD (17.78+0.055) was the lowest. 
Using the Ct values of each experimental sample, 
we draw a box-plot for candidate HKGs (Figure 2). 
The results showed that mean Ct values of HPRT 
(untreated cells: 22.55+ 0.12; RA-treated cells: 
22.51+0.18), ACTB (untreated cells:17.785+0.05; 
RA-treated cells: 17.3340.17) and B2M (untreated 
cells: 25.85+ 0.02; RA-treated cells: 2640.06) did 
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not change significantly after treatment with RA. 
Calculating fold change in expression of the HKGs 
in the treated versus the untreated samples showed 
RA treatment changed the expression of B2M, 
HPRT, GAPDH by 0.86, 1.03 and 1.37 fold. 
Therefore, they can be regarded as reliable 
reference genes for expression analysis of target 
genes in this experiment. On the contrary, 
expression of GAPDH in NCCIT cells 
(Ct:19.36+40.28) reduced significantly (around 
678-fold) after RA treatment (Ct: 28.94+0.18) and 
so this gene is not suitable reference gene for 
expression analysis of target genes after exposure 
to NCCIT and likely any cell lines to RA. 


ACTB 


nue nce 
COS TT 


ONCCIT cells M RA-induced differentiated cells 


GAPDH 


—aa 


ONCCIT cells M RA-induced differentiated cells 


Figure 2. Box-plot based on Ct of 4 candidate HKGs in untreated and RA-treated NCCIT cells. The block diagram 
represents the quartile range (25th to 75th percentiles) of Ct values; the ‘x’ in the box depicts the median of the Ct 
value; underline and overline are determined by the minimum and maximum values of the Ct value. 


This was the first study that analyzed 
expression stability of some HKGs during RA- 
induced differentiation of pluripotent embryonal 
carcinoma NCCIT cells. In another study, stability 
of various HKGs_ during RA-induced 
differentiation of human embryonic stem (ES) cell 
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was examined and £2M, ribosomal protein 
L13A and Alu repeats were found to be the most 
stable for this experimental set-up (Vossaert et al., 
2013). According to another report, both HPRT 
and beta-tubulin mRNA levels varied markedly in 
spontaneously differentiating and growth factor- 
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supplemented (TGF-beta) ES cell cultures, while 
GAPDH expression remained relatively constant 
(Murphy and Polak, 2002). 

RA mediated differentiation and exit from 
the pluripotent state in NCCIT cells were 
confirmed via expression analysis of important 
pluripotency markers including OCT4, NANOG 
and SOX2 using qRT-PCR. The transcription 
factors OCT4, SOX2, and NANOG are important 
factors for maintaining pluripotency and self- 
renewal of pluripotent cells (Chambers and 
Tomlinson, 2009). In many researches, down 
regulation of OCT4, NANOG and SOX2 were 
reported during differentiation of pluripotent cells 
and exit from pluripotency state (Rassouli et al., 
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2013; Soltanian and Dehghani, 2018; Soltanian et 
al., 2014; Soltanian et al., 2020; Stevanovic, 2003). 

Since, B2M, ACTB and HPRT showed stable 
expression levels during RA treatment, they could 
be used as reference genes for expression analysis 
of pluripotency factors between untreated and RA- 
treated NCCIT cells. Results showed that after RA 
treatment, expression of OCT4, NANOG and 
SOX2 decreased to around 79 %, 71 % and 96 % 
of that of untreated NCCIT cells (Figure 3). 
Therefore, according to these results, NCCIT cells 
undergo differentiation by RA treatment. 


RA induced differentiated cells 
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Figure 3. Quantitative RT-PCR analysis of OCT4, NANOG and SOX2 gene expression after treatment of 
NCCIT with RA. NCCIT cells were treated with RA (10 uM) for 14 days and expression of OCT4, NANOG 
and SOX2 were measured and compared to untreated cells as control by real time PCR. The graph represents 
the mean data+SD (error bar) of at least three independent experiments. *p<0.01 against control. 


In conclusion, selection of a suitable HKG is 
essential to normalize real-time PCR data. 
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However, the choice of HKGs should be tailored to 
the nature of the study. For example, although 
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B2M, ACTB, HPRT and GAPDH are widely used 
and accepted as reference genes for analysis of 
gene expression, according to some researches, 
their expression varies in different situation. 
Therefore, expression level stability of HKGs must 
be confirmed in each study (Guénin et al., 2009; 
Zhang et al., 2005). Furthermore, it has been 
strongly proposed that more than one stable 
expressed reference gene should be used to avoid 
misinterpretation of gene expression data 
(Hamalainen et al., 2001; Zhong and Simons, 
1999). 


Acknowledgement 

The authors thank Vice Chancellor for Research 
and Technology, Shahid Bahonar University of 
Kerman, Kerman, Iran for financial support. 


Conflict of Interest 
The authors declare that they have no conflict of 
interest 


References 


Bustin S. A. (2000) Absolute quantification of mRNA 
using real-time reverse transcription polymerase chain 
reaction assays. Journal of molecular endocrinology 
25:169-193. 


Chambers I. and Tomlinson S. R. 
transcriptional foundation of 
Development 136:2311-2322. 


(2009) The 
pluripotency. 


Chen J., Rider D. A. and Ruan R. (2006) Identification 
of valid housekeeping genes and antioxidant enzyme 
gene expression change in the aging rat liver. The 
Journals of Gerontology Series A: Biological Sciences 
and Medical Sciences 61:20-27. 


de Kok J. B., Roelofs R. W., Giesendorf B. A., Pennings 
J. L., Waas E. T., Feuth T., Swinkels D. W. and Span P. 
N. (2005) Normalization of gene expression 
measurements in tumor tissues: comparison of 13 
endogenous control genes. Laboratory investigation 
85:154-159. 


Dheda K., Huggett J. F., Bustin S. A., Johnson M. A., 
Rook G. and Zumla A. (2004) Validation of 
housekeeping genes for normalizing RNA expression in 
real-time PCR. Biotechniques 37:112-119. 


Glare E. M., Divjak M., Bailey M. and Walters E. H. 
(2002) B-Actin and GAPDH housekeeping gene 
expression in asthmatic airways is variable and not 
suitable for normalising mRNA levels. Thorax 57:765- 
770. 


http://jcmr.um.ac.ir 


109 


Soltanian and Sheikhbahaei 


Guénin S., Mauriat M., Pelloux J., Van Wuytswinkel O., 
Bellini C. and Gutierrez L. (2009) Normalization of 
qRT-PCR data: the necessity of adopting a systematic, 
experimental conditions-specific, validation of 
references. Journal of experimental botany 60:487-493. 


Haller F., Kulle B., Schwager S., Gunawan B., von 
Heydebreck A., Siiltmann H. and Fiizesi L. (2004) 
Equivalence test in quantitative reverse transcription 
polymerase chain reaction: confirmation of reference 
genes suitable for normalization. Analytical 
biochemistry 335:1-9. 


Hamalainen H., Tubman J., Vikman S., Kyréla T., 
Ylikoski E., Warrington J. and Lahesmaa R. (2001) 
Identification and validation of endogenous reference 
genes for expression profiling of T helper cell 
differentiation by quantitative real-time RT-PCR. 
Analytical biochemistry 299:63-70. 


Janovick-Guretzky N., Dann H., Carlson D., Murphy 
M., Loor J. J. and Drackley J. K. (2007) Housekeeping 
gene expression in bovine liver is affected by 
physiological state, feed intake, and dietary treatment. 
Journal of Dairy Science 90:2246-2252. 


Jin B., Ernst J., Tiedemann R. L., Xu H., Sureshchandra 
S., Kellis M., and et al. (2012) Linking DNA 
methyltransferases to epigenetic marks and nucleosome 
structure genome-wide in human tumor cells. Cell 
reports 2:1411-1424. 


Kooistra S. M., Thummer R. P. and Eggen B. J. (2009) 
Characterization of human UTF1, a_ chromatin- 
associated protein with repressor activity expressed in 
pluripotent cells. Stem cell research 2:211-218. 


Murphy C. L. and Polak J. M. (2002) Differentiating 
embryonic stem cells: GAPDH, but neither HPRT nor 
B-tubulin is suitable as an internal standard for 
measuring RNA levels. Tissue engineering 8:551-559. 


Nygard A.-B., Jorgensen C. B., Cirera S. and Fredholm 
M. (2007) Selection of reference genes for gene 
expression studies in pig tissues using SYBR green 
qPCR. BMC molecular biology 8:67. 


Park S.-W., Do H.-J., Ha W. T., Han M.-H., Song H., 
Uhm S.-J., and et al. (2014) Differential expression of 
ETS family transcription factors in NCCIT human 
embryonic carcinoma cells upon retinoic acid-induced 
differentiation. Biological and Pharmaceutical Bulletin 
37:659-665. 


Radoni¢ A., Thulke S., Mackay I. M., Landt O., Siegert 
W. and Nitsche A. (2004) Guideline to reference gene 
selection for quantitative real-time PCR. Biochemical 
and biophysical research communications 313:856-862. 


Rassouli F. B., Matin M. M., Bahrami A. R., 
Ghaffarzadegan K., Cheshomi H., Lari S., Memar B. 


Journal of Cell and Molecular Research (2021) 12 (2), 104-110 


and Kan M. S. (2013) Evaluating stem and cancerous 
biomarkers in CD15+ CD44+ KYSE30 cells. Tumor 
Biology 34:2909-2920. 


Richly E., Chinnery P. F. and Leister D. (2003) 
Evolutionary diversification of mitochondrial 
proteomes: implications for human disease. Trends in 
Genetics 19:356-362. 


Selvey S., Thompson E. W., Matthaei K., Lea R. A., 
Irving M. and Griffiths L. R. (2001) B-Actin—an 
unsuitable internal control for RT-PCR. Molecular and 
cellular probes 15:307-311. 


Soltanian S. and Dehghani H. (2018) BORIS: a key 
regulator of cancer stemness. Cancer cell international 
18:1-13. 


Soltanian S., Dehghani H., Matin M. M. and Bahrami A. 
R. (2014) Expression analysis of BORIS during 
pluripotent, differentiated, cancerous, and _ non- 
cancerous cell states. Acta Biochim Biophys Sin 
46:647-658. 


Soltanian S., Sheikhbahaei M. and Ziasistani M. (2020) 
Phytol Down-Regulates Expression of Some Cancer 
Stem Cell Markers and Decreases Side Population 
Proportion in Human Embryonic Carcinoma NCCIT 
Cells. Nutrition and Cancer: 1-14. 


Stevanovic M. (2003) Modulation of SOX2 and SOX3 
gene expression during differentiation of human 
neuronal precursor cell line NTERA2. Molecular 
biology reports 30:127-132. 


Turabelidze A., Guo S. and DiPietro L. A. (2010) 
Importance of housekeeping gene selection for accurate 
reverse transcription-quantitative polymerase chain 
reaction in a wound healing model. Wound repair and 
regeneration 18:460-466. 


Vandesompele J., De Preter K., Pattyn F., Poppe B., Van 
Roy N., De Paepe A. and et al. (2002) Accurate 
normalization of real-time quantitative RT-PCR data by 
geometric averaging of multiple internal control genes. 
Genome biology 3:1-12. 


Vossaert L., O’Leary T., Van Neste C., Heindryckx B., 
Vandesompele J., De Sutter P. and et al. (2013) 
Reference loci for RT-qPCR analysis of differentiating 
human embryonic stem cells. BMC molecular biology 
14:21. 


Wu Y.-Y. and Rees J. L. (2000) Variation in epidermal 
housekeeping gene expression in different pathological 
states. Acta dermato-venereologica 80: 2-3 


Zhang X., Ding L. and Sandford A. J. (2005) Selection 
of reference genes for gene expression studies in human 
neutrophils by real-time PCR. BMC molecular biology 
6:4. 


http://jcmr.um.ac.ir 


Soltanian and Sheikhbahaei 


Zhong H. and Simons J. W. (1999) Direct comparison 
of GAPDH, B-actin, cyclophilin, and 28S rRNA as 
internal standards for quantifying RNA levels under 
hypoxia. Biochemical and _ biophysical research 
communications 259:523-526. 


Open Access Statement: 

This is an open access article distributed under the 
Creative Commons Attribution License (CC-BY), 
which permits unrestricted use, distribution, and 
reproduction in any medium, provided the original 
work is properly cited. 


110 


Journal of Cell and Molecular Research (2021) 12 (2), 104-110 (S1) Soltanian and Sheikhbahaei 


Supplementary Materials: 


A 
Untreated NCCIT cells RA-induced differentiated NCCIT cells 
B 
Sample Concentration 260/280 260/230 
Undifferentiated 1000 pg/ml 2.1 18 
NCCIT cells 
RA- induced 1000 pg/ml 2.2 1.9 


differentiated cells 


Figure S1. Agarose gel electrophoresis of total RNA isolated from undifferentiated and RA- induced 
differentiated NCCIT cells. A) Intact total RNAs have sharp 28S and 18S rRNA bands. B) Nanodrop analysis 
of RNA samples isolated from undifferentiated and RA- induced differentiated NCCIT cells. 
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Figure $2. Melt-curve analysis of products from SYBR Green assay for each primer sets. 
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Table S1. The man Ct value of three biological repeats 
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Gene Name B2M HPRT GAPDH ACTB 
Ct values 25.88 22.69 19.58 17.84 
(NCCIT 25.83 22.5 19.04 17.73 
untreated 25.85 22.46 19.48 17.785 
celCs) 

Ct values 25.99 22.69 29.12 17.16 
(RA- 26.72 22.33 28.76 17.51 
treated 26.09 22.57 28.94 17.335 
NCCIT 

celCs) 
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